External field control of donor electron exchange at the Si/SiC>2 interface 
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We analyze several important issues for the single- and two-qubit operations in Si quantum com- 
puter architectures involving P donors close to a SiC>2 interface. For a single donor, we investigate 
the donor-bound electron manipulation (i.e. 1-qubit operation) between the donor and the inter- 
face by electric and magnetic fields. We establish conditions to keep a donor-bound state at the 
interface in the absence of local surface gates, and estimate the maximum planar density of donors 
allowed to avoid the formation of a 2-dimensional electron gas at the interface. We also calculate 
the times involved in single electron shuttling between the donor and the interface. For a donor 
pair, we find that under certain conditions the exchange coupling (i.e. 2-qubit operation) between 
the respective electron pair at the interface may be of the same order of magnitude as the coupling 
in GaAs-based two-electron double quantum dots where coherent spin manipulation and control 
has been recently demonstrated (for example for donors ~ 10 nm below the interface and ~ 40 nm 
apart, J ~ 10 -4 meV), opening the perspective for similar experiments to be performed in Si. 

PACS numbers: 03.67.Lx, 85.30.-z, 73.20.Hb, 85.35.Gv, 71.55.Cn 
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I. INTRODUCTION 

Doped Si is a promising candidate for quantum com- 
putingi due to its scalability properties, long spin coher- 
ence times £&d£&il and the astonishing progress on Si 
technology and miniaturization in the last few decades 
(Moore's law). The experimental production of a work- 
ing qubit depends on precise positioning (of the order 
of A} 8 ' 9 of donors in Si and the quantum control of the 
donor electrons by local gates placed over an oxide layer 
above the donors. The required accuracy in donor po- 
sitioning has not been yet achieved, although there are 
increasing efforts in this direction, using top-down tech- 
niques, i.e. single ion implantation (with tens of nm accu- 



racy) 



10.11.12 



and bottom up techniques, i.e. positioning 



of P donors on a mono-hydride surface via STM (with 
1 nm accuracy) with subsequent Si overgrowth i 13 : 14 

In the original doped Si based quantum computer pro- 
posal, 1 illustrated in Fig. [I] the qubits are the donor nu- 
clear spins, and the hyperfine interaction between these 
and the donor electron spins is used to perform single 
qubit operations (rotations). The strength of the hyper- 
fine interaction is manipulated by local surface gates, the 
so called A-gates, which move the electron between the 
donor and an interface with SiC>2. Exchange between 
neighboring donors, tuned by surface 'exchange' gates 
(J-gates), would control two-qubit operations. Exchange 
gates were originally proposed for a double quantum dot 
geometry in GaAs. 15 Related proposals in Si use the elec- 
tron spin as qubi t 16 ' 17 ! 18 or the electron charged Charge 
coherence in Si is much shorter (~ 200 ns)2£i than the spin 
coherence Ti ~ 1 ms, which can be further enhanced by 
isotopic purification^**"! making spin qubits in general 
more attractive than charge qubits for actual implemen- 
tations. On the other hand, direct detection of a sin- 
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FIG. 1: (Color online) Schematic view of Kane's quantum 
computer. Donors D are a distance d from the Si/Si02 (001) 
interface and a distance R from each other. Surface A- and 
J-gates control one and two qubit operations. In the present 
study, we consider uniform electric F and magnetic B fields 
applied in the z-direction. 



gle spin is a very difficult task ) 21 : 22 : 23 while a fraction 
of a single electron charge can be easily detected with 
state-of-the-art single electron tunneling (SET) devices. 
As a result, ingenious spin-to-charge conversion mech- 
anisms, that would allow the electron spin state to be 
inferred according to the absence or presence of charge 
detected by an SET at the surface, have been proposed, 
e.g. Refs. [IHHlil 

Doped Si has two main advantages over GaAs quan- 
tum dots: (i) The much longer spin coherence times, that 
can be enhanced by isotopic purification (note that all 
isotopes of Ga and As have nuclear spins so the spin co- 
herence time in GaAs cannot be improved via isotopic 
purification), and (ii) the identical Coulomb potentials 
created by donors as opposed to variable quantum dot 
well shapes produced by surface gates on a 2-dimensional 
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electron gas (2DEG). Despite this latent superiority of 
Si, progress in GaAs has been much faster ) 23 i 28 i 29 ' 30 i 31 
in particular due to the fact that the electrons, being 
at the device surface, are easier to manipulate and de- 
tect. Another Si handicap is that the exchange between 
donors in bulk Si oscillates, changing by orders of mag- 
nitude when the relative position of neighboring donors 
changes by small distances ~ This is caused by in- 

terference effects between the six degenerate minima in 
the Si conduction band. However, as discussed below, 
this degeneracy is partially lifted at the interface, thus 
the oscillatory behavior may not represent such a severe 
limitation for interface states as compared to the donor 
bulk states. 

In the following, we analyze the manipulation of donor 
electrons close to a Si/Si02 interface by means of exter- 
nal uniform electric and magnetic fieldsj 33 ' 34 In Sec. |TT] 
we introduce the model for an isolated donor and discuss 
the interface and the donor ground states, which are cal- 
culated variationally. We also analyze the shuttling be- 
tween the interface and the donor, including the effect of 
a magnetic field. In Sec. IIHI we study the conditions to 
avoid the formation of a 2DEG at the interface, and dis- 
cuss the advantages and actual feasibility of performing 
two-qubit operations at the interface. A summary and 
conclusions are given in Sec. HVl 



II. SINGLE DONOR 
A. Model 
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FIG. 2: (Color online) Double-well potential formed by the 
donor nuclear charge (Wd), the applied electric field, and the 
interface with the oxide (considered to be oo, i.e., impene- 
trable). The interface well Wj consists of the electric field 
and the interface in the ^-direction and also includes the nu- 
clear charge and its image in the xy-plane [(see Fig. [3ja)]. 
The dashed lines represent the decoupled ground eigenstates 
in each well ipz and %po- The thin horizontal lines indicate 
the expectation value of the energy in each well. The po- 
tential, wave-functions and energies depicted correspond to 
d = 30 nm and F = 13.5 kV/cm w F c (30nm). 



We consider initially a single donor a distance d from 
a Si/Si02 (001) interface. As a simple model for the 
A-gate effects, a uniform electric field is applied in the 
z-direction, as illustrated in Fig. [TJ The inter-donor dis- 
tance R is assumed to be large enough so that each 
donor can be treated as an isolated system. The con- 
duction band of Si has six equivalent minima located 
along the A lines. As discussed below, it is reason- 
able to treat this system within the single- valley effective 
mass approximatio n 35 i 36 i 37 i 38 i 39 leading to the following 
Hamiltonian for the donor electron 



H = T + eFz- 



e 2 Q 



e 2 Q 



esir e Siv /p 2 + {z + 2d) 2 4e S i(z + d) 

We also consider an applied magnetic field along z: 
In this case the vector potential A = B(y,—x,0)/2 
is included in the kinetic energy term, T = 
^ v =x, y ,z K 2 /^) [id/dv + eA v /(hc)f. The effective 
masses in Si are m x = m y = m± = 0.191m, and m z = 
rn 1 1 = 0.916 to. The second term is the electric field linear 
potential, the third is the donor Coulomb potential, and 



the last two terms (with p 



+ y ) take account of 



the charge images of the donor and the electron, respec- 
tively. Q = (e S i0 2 - esi)/(esi0 2 + ^Si), where e S i = 11.4 
and esio 2 = 3.8. In this case Q < and, therefore, the 



images have the same sign as the originating charges. In 
rescaled atomic units, a* = Ti 2 esi/ m _\_e 2 = 3.157 nm and 
Ry* = m^e 4 /2h 2 e 2 Si = 19.98 meV, and the Hamiltonian 
is written 



H = - 



+ neFz - 



d 2 d 2 

dx 2 dy 2 
2 



d 2 

oz z 



1 2 2 



2Q 



i/J-(yd x - xdy) 
Q 



r y/p? + (z + 2d) 2 2{z + d) ' 



(2) 



where 7 = m±/m^, fi 2 = a* 4 /A|j with Xb — yh/eB the 
magnetic length, k = 3.89 x 10~ 7 e| i (to/tojJ 2 cm/kV, 
and the electric field F is given in kV/cm. 

The system under study consists of a shallow donor, 
P in particular, immersed in Si a distance d from a Si02 
barrier, which we assume to be infinite (impenetrable). 
When no external field is present, the electron is bound to 
the donor potential well Wd ■ When an electric field F is 
applied in the z-direction, a triangular well Wi is formed 
next to the interface. The interface well Wi also includes 
the donor and its image Coulomb potentials at the in- 
terface (z — —d) which, under special circumstances dis- 
cussed below, are strong enough to localize the electron 
in the xy-plane. The interface well Wi and the donor 
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Coulomb potential Wd form an asymmetric double well, 
as illustrated in Fig. [2] 

The Hamiltonian in Eq. j2]) is solved in the basis 
formed by ipo and ipi, which are the ground eigenstates 
of each of the decoupled wells Wd and Wi . The Hamil- 
tonian is written 



H = T + Vi + V D 



(3) 



where the last term is added to avoid double counting of 
the impurity Coulomb potential at the interface included 
both in Vi and Vd, which are defined as 



Vi 



, Q-i 

' VP 2 + d2 



KeFz - 



Q 



2(z + d) 



and 



Vd = — 



(4) 



(5) 



The first term in Eq. ^ is the sum of the donor Coulomb 
potential and its image charge potential at the interface. 

The Hamiltonian in the non-orthogonal basis {ipD, 4>i} 
reads 



Hdd Hid 
Hid Hjj 



(6) 



where Hab = (^a\H\iPb) with A,B = I,D. Diago- 
nalization gives the two eigenstates ty + and with 
eigenenergies E + and E~ which show anticrossing behav- 
ior with a minimum gap when Hdd — Hn. This point 
defines the characteristic field F c , illustrated in Fig. [21 
which is relevant for the tunneling process discussed in 
detail in Sec. iHDl 



B. Interface state tpi 

It is convenient to write the interface potential Vj as a 
sum of purely z- and purely p-dependent terms: 



Vi 

Vj" 



Vf 



V? + Vf 
neFz 



2(z + d) 



1 



vV + d 2 



(7) 
(8) 

(9) 



The Vf component is the triangular well plus the elec- 
tron image charge potential, while the Vf component is 
the sum of the impurity and its image potential at the 
interface. Vf is plotted in Fig. [3£a) for three different 
values of d. Curves corresponding to the parabolic ap- 
proximation of the potential, 
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FIG. 3: (Color online) (a) Confining potential at the inter- 
face xy-plane due to the Coulomb attraction of the donor 
nucleus and its image Vf = 2(Q - l)/vV + d 2 - Dotted lines 
correspond to the harmonic approximation of the potential 
Vparab(p) = {Q - l)(2/rf - p 2 / 'd 3 ) . (b) Ground state energy 
E p (solid lines) and first excited state energy E' p at the in- 
terface (dashed lines). The thick lines correspond to the po- 
tential Vf while the thin lines correspond to the parabolic 
approximation. As d increases, the solution for the parabolic 
approximation approaches the one for the Vf potential, as 
expected. 



are also shown, and it is clear that the harmonic approx- 
imation works better for the larger distances d. Elec- 
tron confinement at the interface xy-planc is provided by 
Vf, while the uniform electric field and the oxide confine 
the electron along the z-direction. For certain values of 
d and R, Vf is deep enough to localize the individual 
donor electrons (with no need of local A-gates) and keep 
them from forming a 2DEG at the interface: The neces- 
sary conditions are discussed in Sec. IIII Al The spacial 
localization of the electrons at the interface is a neces- 
sary condition for Si-based quantum computing if qubit 
read-out takes place at the interfaced 

The ground state at the interface is calculated by solv- 
ing Hi = T+Vi variationally with a separable trial func- 
tion 



4>i = f(z)g(p) 



(ii) 



For the z-part we use 



(2£+l)/2 

f( z ) = (z + dy e - a( - Z+d)/2 (12) 



for z > —d. The infinite barrier at the interface is taken 
into account by forcing the ground state to be zero at 
the interface, so f(z) = for z < —d. a is a variational 
parameter that minimizes the contribution to the energy 
(f(z)\T z + Vf\f(z)) . The most suitable value for I is cho- 
sen by comparing E z = (f(z)\T z + neFz\f(z)) and f(z) 
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FIG. 4: (Color online) (a) E z for different trial functions f(z) 
compared with the exact solution of the infinite triangular 
potential.— (b) Resulting wave-function for F = 50 kV/cm. 
z' — z + d. Both the energy and wave-function for I — 2 give 
the best match to the exact solution. 



with the exact solution of an infinite triangular well 4 ^ 

/ oxact (z') = Ai 



i2m 7 eF 



eF 



(13) 



where Ai is the Airy function, z' = z + d, and Eq is the 
ground state energy 



Eq — E c 



' (0.7587 neF) 2 . (14) 



2m 



The results are shown in Fig. [4] for I = 1, 2, and 3. 
I = 2 gives the best agreement with the exact solution for 
both the energy and the wave-function and is therefore 
adopted in what follows. 

For the p-part we use the ansatz 



9 (P) = -^e-^/i 



(15) 



with the variational parameter (3 calculated by minimiz- 
ing E p = (g{p)\TP + Vf\g{p)). We have checked that 
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FIG. 5: (Color online) Typical confinement lengths of the 
interface state in the xy-plane 1//3 and z-direction 1/a. a 
depends on the value of the electric field applied: a c shown 
here corresponds to F c (d). 



this gaussian form gives lower energy than an exponen- 
tial e _I,p / 2 (as a reminiscent of the donor wave-function) 
for distances d > 1 nm. 

In Fig. ODJb) we plot the energy E p for the ground state 
and the first excited state (g'(p) oc xgpi{p)) for both the 
variational solution adopted here and the parabolic ap- 
proximation of Vf . The parabolic approximation gives 
an underestimation of the binding energies and diverges 
at short distances d (not shown). To guarantee that the 
electron remains bound, and at the ground state, the 
operating temperature has to be lower than the energy 
difference between the ground and the first excited states 



k B T <C min(\E p \ 



Ep). For d = 30 nm, the exci- 



tation gap Ep — Ep ~ 1 meV. This limits the operat- 
ing temperature to a few Kelvin (in current experiments, 
temperatures as low as 0.1 K are being used). 41 ! 42 

The inverse of the variational parameters, 1/a and 
1//3, are proportional to the confinement lengths in the 
z-direction and the xy-plane respectively. Both depend 
on the distance d of the donor from the interface as shown 
in Fig. [5j and 1/a also depends on the value of the ap- 
plied electric field, being larger for smaller F. Fig. Ogives 
the value of 1/a for the electric field F c (d) at which the 
expectation values of the energies of the states at the in- 
terface and at the donor are degenerate. The expectation 
value of the position of the electron in the z-direction is 
J_ d {z + d)f 2 (z) dz = 5/a from the interface. This value 
is small compared to d (for d = 30 nm, 5/a = 5 nm), 
justifying the validity of the two-well approach we are 
using to solve the Hamiltonian. 

At the interface, the z-valleys' energy is lower than the 
xy- valleys'. This is straightforward to show for an infi- 
nite triangular potential 4 ^ in which E z (m±) / E z {m\\) = 
1 / 3 = 1.68637 [see Eq. dHJ)] . The difference be- 
tween the levels depends on the electric field as F 2 / 3 . 
AE = E z (m±^) — E z (m\\) is shown in Fig. [7] (a) (right 
axis). For a field F = 5 kV/cm, which is small for our 
interests, the splitting is Ai? = 2.5 meV which corre- 
sponds to T ~ 30 K. If a magnetic field is applied, the 
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z-levels increase their energy faster than the ly-levels un- 
til they cross. However, this crossing happens at a very 
large magnetic field B > 20T (shown in the left axis of 
Fig. [7] (a) as a function of the electric field F). Therefore, 
for the range of parameters of interest here, the ^-valleys 
are always the ground state at the interface. We point 
out that it is experimentally established that, in MOS- 
FET geometries equivalent to the one studied here, the 
interface ground state is non-degenerate, with a 0.1 meV 
gap from the first excited state.— This is well above the 
operation temperatures in the quantum control experi- 
ments investigated in the present context. 

The magnetic field has two main effects on the states: 
(i) the electron gets more confined in the direction par- 
allel to the interface and, consequently, (ii) its kinetic 
energy increases. The effect of the magnetic field is 
strongest for the less confined wave-functions, which cor- 
respond to the larger cf s. We can quantify the strength 
of this effect by calculating the magnetic field B c that is 
needed to get a magnetic length As of the same order of 
the confinement length in the plane parallel to the inter- 
face: for a donor a distance d = 30 nm, 1/0 = 18.5 nm 
and B c ~ 2 T while for d = 15.8 nm, 1/0 = 12 nm and 
B c ~ 4.5 T. The confinement effect is illustrated in Fig. [6] 
where 1/ 0b as a function of magnetic field for two differ- 
ent values of d is shown. The thick lines correspond to 
the variational solution of minimizing 



dx 2 



dy 2 



1 2 2 

V 1 p 



2Q 



d 2 



(16) 



with trial function ~ exp(— 0gp /2). Closer donors pro- 
duce a larger confinement of the interface electron wave- 
function but the effect of the magnetic field is much more 
dramatic for the donors further away from the interface: 
for d — 30 nm, a magnetic field of 10 T decreases the 
wave- function radius by a 40%. Within the parabolic ap- 
proximation for the interface potential V r para b(p) (which 
overestimates the wave-function confinement) the depen- 
dence of 0b on the field is given by^ 



B = (/5 4 + 



1/4 



(17) 



with values as shown in Fig. [5] (thin lines). 

The increase in energy with magnetic field, as cal- 
culated variationally, is shown in Fig. [7] (b). We ob- 
serve again that the effect of the magnetic field is much 
stronger for the larger values of d. The much smaller 
shift in the donor ground state energy (discussed in the 
next subsection) is also shown for comparison. 




2 4 6 8 10 
Magnetic Field (T) 

FIG. 6: (Color online) 1/0 versus magnetic field for two differ- 
ent interface-donor distances. Thick lines represent the vari- 
ational results and the thin lines correspond to the harmonic 
approximation in Eq. (|17[) . 
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FIG. 7: (Color online) (a) (Right axis) Splitting between the 
z- and the xy-\evels at the interface due to the anisotropic 
mass. (Left axis) Magnetic field at which the xy-levels become 
the ground state at the interface, (b) Energy shift of the 
interface (z-valleys) and donor levels (both z- and xy- valleys) 
under an applied magnetic field. Due to the mass anisotropy, 
the shift with magnetic field for the z-valleys levels is larger 
than for the xt/-levels. At the donor, this implies that the 
degeneracy of the valleys is broken. At the interface, the z- 
valleys are the ground state up to very large values of the 
magnetic field as shown in (a). 



C. Donor state ifin 

The potential Vd consists of the isolated impurity 
Coulomb potential 



The solution of Ho = T + Vd is taken to be of the 
form of the anisotropic envelope wave- function^ multi- 
plied by (z + d) to satisfy the boundary condition at the 
interface ipD\z=-d = 



ip D = N(z + d)e-v /p2/a2+z2/b \ 



( 18 ) where 1/N 2 



•■b(d 2 



(19) 



\be~ 2d / b {\d + b)), for 




d (nm) 

FIG. 8: (Color online) Variational donor ground state energy 
(iI>d\Hd\'>Pd) versus d obtained from our trial function tpo and 
from the trial function used by MacMillen and Landman 38 
(data taken from Table III in Ref. HH). The short line at large 
d represents the value of the ground state energy of an isolated 
donor in bulk Si (within the single valley approximation) ~ 
—31.2 meV.— A "perfectly imaging" plane is used (e xide = 
0), as assumed in Ref. |38| . For the intermediate and large 
values of d of interest here, results from the two approaches 
are in reasonable agreement. 
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Electric field (kV/cm) 

FIG. 9: (Color online) Eigenenergies E + and E~ of the 
Hamiltonian in Eq. © for d — 11 nm. They show an anti- 
crossing behavior with a minimum gap g m i n at F c . Tunneling 
times are related to grain as r = fr/g m in- F c can be determined 
as the field at which Hji = Hdd- 



z > —d. tpD = for z < —d. For 
d >> a,b, ipo reduces to the bulk limit ^ ulk = 

l/(7ra 2 6) exp (^—\/p 2 /a 2 + z 2 /b 2 ^j. a and b are varia- 
tional parameters chosen to minimize the ground state 
energy. Except for the smallest distances d < 2a* ~ 
6 nm, not relevant here, a and b coincide with the Kohn- 
Luttinger variational Bohr radii of the isolated impurity 
(d — > oo) a — 2.365 nm and b = 1.36 nm. 

In Fig. [5] we show the variational results for the ground 
state energy obtained from our trial function tpD- For 
comparison, we also give results obtained through the 
trial function proposed by MacMillen and Landman, 38 
where, aiming at a good description for donors at very 
short distances from the interface (typically smaller than 
the effective Bohr radius a*), a much larger set of vari- 
ational states was used for the expansion of the donor 
state. For this comparison, our results in Fig. [8] corre- 
spond to a "perfectly imaging" plane {e xide = 0), as 
assumed in Ref. [H. The energy depends strongly on d 
for the smaller values of d, and tends to the bulk value 
at long distances. For the intermediate and large values 
of d of interest here, the two approaches are essentially 
equivalent. 

We find that the effect of the external fields on the 
donor state is negligible. For instance, for the largest 
electric fields of interest here (F c at short distances 
d ~ 6 nm) , the energy corresponding to the electric field 
potential is (ip d\kcF 'z\tJj d) = 0.18 Ry*, to compare with 
1.6 Ry* for the isolated donor ground state in bulk. The 
effect of a magnetic field on the electron wave-function 
at the donor is also very small: the field required to 
get a magnetic length of the order of the Bohr radius 
a = 2.365 nm is B c ~ 120 T! The donor ground state 
energy shift due to the magnetic field can be estimated 



byii 

AE = ^Ry* , (20) 

where lengths are given in units of a*, (f 2 ) = 2a 2 /a* 2 
for the z-envelopes, and (r 2 ) = (a 2 + b 2 )/a* 2 for the xy- 
envelopes. The results are shown in Fig. [7]^b) and are 
comparable to the values for the Is orbital of shallow 
donors calculated numerically in Ref. |45j. Note that the 
magnetic field partially breaks the six-valley degeneracy 
due to the different confinement radius of the electron 
wave- function in each of the different valleys. 



D. Shuttling between interface and donor states 

We model the donor electron ionization under an ap- 
plied electric field along z by considering the tunneling 
process from the donor well into the triangular well at the 
Si/SiC"2 interface (see Fig. [2]). The required value of the 
field for ionization to take place may be estimated from 
Ha = Hdd [see Eq. ©]. We call F c the characteristic 
field for which this condition is fulfilled, which is equiva- 
lent to require that the gap between the two eigenenergies 
E + and E~ of H is minimum, as illustrated in Fig. [5] 

Our results for F c versus d are shown by the solid dots 
[labelled (1)] in Fig. [TQl In this figure we also test the 
robustness of our approach, namely using tJjd as given in 
Eq. (fTO|) , by comparing the values of F c obtained assum- 
ing different forms for the donor trial function. Curves 
(2) and (3) correspond to isotropic and anisotropic wave- 
functions respectively, with the same ground state energy 
for the electron at the donor as obtained from in 
Eq. (HHJ), ~ —31.2 meV. Note that they compare very well 
with curve (1). Curve (4) corresponds to a tight-binding 
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FIG. 10: (Color online) Comparison of F c versus d for the 
single donor problem (P) obtained from different forms for 
ij)D, Label (1) corresponds to ipo with anisotropic mass as 
defined in Eq. (|19|l . Label (2) assumes isotropic effective mass 
m* = 0.29819m, chosen so that the ground state energy is 
—31.2 meV, the same as for the single valley approximation. 

Label (3) corresponds to ipt> ~ exp ^— p 2 ja 2 + z 2 /b 2 ^j , 

and the same ground state ener gy. Label (4) reproduces 
the tight-binding results in Ref. l4q . where the six valley 
degeneracy is considered, leading to the ground state energy 
of —45 meV. The latter coincides with the experimental 
binding energy for a P donor in Si. Label (5) considers 
tpo ~ exp(— r/a) with isotropic mass m* = 0.43m so that 
the ground state energy is —45 meV. 



result^ in which the six-valley degeneracy of Si is incor- 
porated. Although all curves arc qualitatively similar, 
curve (4) is shifted towards larger fields. The origin of 
this shift is investigated by considering an isotropic trial 
function whose parameters have been chosen to give a 
ground state energy ~ —45 meV, and we note that the 
results, given in curve (5), compare very well with those 
in curve (4). We conclude that the shift in the value of F c 
when the six-fold degeneracy of the Si conduction band 
is considered is mainly due to the fact that the ground 
state energy at the donor in the single valley approxima- 
tion [curve (1)] is ~ —31.2 meV while it is ~ —45 meV 
when the intervalley coupling is included [curve (4)]. In 
the following we use ipo as defined in Eq. (|19p , keeping in 
mind that electric field values are bound to be somewhat 
underestimated. 

We may picture the electron shuttling between the two 
wells under applied electric and magnetic fields by calcu- 
lating the expectation value of its position along z at 
the ground state, (z) = (\I> - [zl^ - ). The results for 
d = 20.5 nm are shown in Fig. [IT] , where the horizon- 
tal lines mark the position of the interface. The distance 
between (z) and the interface tends to 5/ a for F > F c , 
where a also depends on F. In Fig. Illf a) we show pic- 
torially how the electron would evolve from the donor 
to the interface well when an electric field is applied. 



20 30 1 2 3 4 
Electric field (kV/cm) Magnetic field (T) 

FIG. 11: (Color online) Expectation value of the electron 
z-coordinate (z) = for d = 20.5 nm versus (a) 

electric field and (b) magnetic field [F — F c + 60 V/cm]. The 
horizontal lines represent the relative position of the interface. 
The donor is at (z) = 0. The electric field moves the electron 
from the donor to the interface while a parallel magnetic field 
takes the electron back to the donor. 



At small values of F the electron is at the donor well, 
(z) ~ and w ipr>- The center of mass is slightly 
shifted from the donor site due to the factor (z + d) in 
ipD. Above F c , the electron eventually tunnels to the 
interface. Starting with the electron at the interface in 
a near-degeneracy configuration (F > F c ), a relatively 
modest magnetic field can cause the electron to move in 
a direction parallel to the field and against the electric 
field, as shown in Fig. [TTT b). This is due to the much 
larger shift of the interface state energy with magnetic 
field compared to the shift of the donor ground state en- 
ergy [see Fig.[7fb)]. This behavior characterizes electrons 
originating from the donors, and not other charges that 
the SET may detect, like charges in metallic grains on 
the device surfaced The combination of parallel electric 
and magnetic fields constitutes therefore a valuable ex- 
perimental setup to investigate whether charge detected 
at the interface actually originates from a donor^i 

A key parameter determining the feasibility of quan- 
tum computation in the doped Si architecture is the time 
required to shuttle the electron between the donor and 
the interface. This time should be orders of magnitude 
smaller than the coherence time to allow for many op- 
erations and error correction while coherent evolution of 
the qubit takes place. The tunneling process conserves 
the spin, but coherence would be lost for orbital/charge 
degrees of freedom. Therefore, if quantum information is 
stored in a charge qubit, the electron should evolve adia- 
batically from the donor to the interface, while tunneling 
would be acceptable for spin qubits. In an adiabatic pro- 
cess the modification of the Hamiltonian (for instance, 
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FIG. 12: (Color online) Donor ionization tunneling (solid line) 
and adiabatic passage (dashed line) times versus d. 



when an external field is applied) is slow enough that the 
system is always in a known energy eigenstate, going con- 
tinuously from the initial to the final eigenstate— Here, 
we calculate both tunneling and the adiabatic passage 
times. 

We estimate the tunneling time from the minimum gap 
5min between the two eigenvalues E + and E~ (see Fig.O 
via the uncertainty relation r = h/g m i n . The adiabatic 
time is calculated a o 46 i 48 r a = h\e\F max d/g^ nin and is or- 
ders of magnitude larger than the tunneling time. -F m ax 
is chosen so that the electron is at the interface 'F - w ipi- 
The results for the tunneling and adiabatic passage times 
(for F max = 2F C ) are shown in Fig.[T2] The times depend 
exponentially on the distance d. Tunneling times range 
from 0.1 ps for d — 6 nm to 10 ns for d = 38 nm. Adi- 
abatic times range from 1 ps for d = 6 nm to 100 ns for 
d = 20 nm and get very large at longer distances. These 
times are to be compared to the experimental values of 
spin coherence and charge coherence respectively. 

Spin dephasing in Si is mainly due to dipolar fluctua- 
tions in the nuclear spins in the system, which produce 
a random magnetic field at the donor electron spin. The 
spin dephasing time in bulk natural Si is Ti ~ 1 ms^i 4 ^ 
Natural Si is mostly composed of 28 Si (no nuclear spin), 
with a small fraction (4.67%) of 29 Si (nuclear spin 1/2), 
therefore, X2 can be dramatically improved through iso- 
tope purification 2 -! 4 .^ up to 100 ms or longer in bulk. 
Moreover, it has been recently proposed that the spin 
dephasing times can be arbitrarily prolonged by apply- 
ing designed electron spin resonance pulse sequences^ On 
the other hand, the closeness of a surface or interface can 
reduce the spin dephasing times4^ The tunneling time 
calculated here is orders of magnitude smaller than the 
bulk T2 ~ 1 ms for natural Si: In the worst case scenario 
(long distances d ~ 40 nm) T2/T ~ 100. Therefore, we 
expect tunneling times to be always much shorter than 
decoherence times even in the presence of an interface. 

Charge coherence time is considerably shorter than 
spin coherence time, since charge couples very strongly 
to the environment through the long range Coulomb in- 
teraction. The main channels of decoherence are charge 
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FIG. 13: (Color online) The main panel gives the maximum 
donor planar density n max for which electrons drawn to the 
interface remain localized around the donor region, and do 
not form a 2DEG. It is assumed that all donors are at the 
same distance from the interface d. n max is estimated from 
the criterion that the barrier between neighboring wells at 
the interface is above Edw given in Eq. (1 2 2 [I . Results are 
shown for B = T and B = 10 T. The inset shows the 
double well potential parallel to the interface created by two 
donors located a distance d — 6.3 nm from the interface and 
separated by R = 28 nm. The expectation value of the energy 
Bdw (given by the dashed lines for B = 0) is also lower than 
the single-donor well ground state energy. 



fluctuations and electron-phonon interactionsi 50 ' 51 The 
charge coherence time has been measured to be ~ 200 ns 
for Si quantum dots surrounded by oxide layers.— This 
number has to be compared to the calculated adiabatic 
times, which are much longer than the tunneling times. 
Therefore, for charge qubits to be realizable in the con- 
figuration discussed here, the donor-interface distance d 
has to be limited to a maximum of 20 nm. 



III. DONOR PAIR 

A. Planar density 

We estimate the maximum planar density of donors 
("■max) allowed to avoid the formation of a 2DEG as 

"max 

= l/i? 2 lin , where i? m in is a minimum distance be- 
tween two donors which is calculated as follows: We as- 
sume two donors located at the same distance d from the 
interface and a distance R apart. The resulting double 
well potential along the interface xy-plane is 



Vtw{x,v) 



2(Q-1) 



y/(x + R/2) 2 +y 2 +di 

2(Q-1) 
y/{x - R/2) 2 + y 2 + d 2 



(21) 
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as illustrated in the inset of Fig. [13] for d = 6.3 nm and 
R = 28 nm. We adopt two different criteria to estimate 
-Rmin- (i) The first one requires i? m i n = 2/(3 where 2/(3 is 
the width of the gaussian g(p) (see Fig. [SJj. For instance, 
for d — 30 nm, this gives i? m in ~ 40 nm leading to n max ~ 
6 x 10 10 cm~ 2 , while for d = 10 nm, i? m i n ~ 18 nm and 
"■max ~ 3 x 10 11 cm -2 , (ii) The second criterion, which 
we find to be slightly more restrictive, requires a high 
enough barrier within the double well, and is given by 

Ebw = (ikl Huw hM < V DW (0, 0) , (22) 

where i? m i n corresponds to the equality condition. -S/dw 
is the double-well 2-dimensional Hamiltonian 

H DW =T x +T y + V DW (x,y), (23) 

with T x and T y the kinetic energy terms. Vdw(0,0) is 
the maximum height of the inter-well barrier, which is re- 
quired to be above the single-particle expectation value 
of the energy -E"dw- The maximum planar density es- 
timated from Eq. (|22p is shown in the main panel of 
Fig. E2 For instance, n max (d = 30 nm) « 10 10 cm~ 2 
obtained from i? m i n w 88 nm. n max is larger for the 
donors closer to the interface. For instance, n max (d — 
10 nm) w 6 x 10 10 cm" 2 (i? min w 38 nm). 

As shown in Fig.[BJ 1/ [3b decreases with a perpendicu- 
lar magnetic field, hence increasing the maximum planar 
density. The first criterion for the maximum planar den- 
sity gives, for d = 30 nm and B — 10 T, i? m i n s» 22 nm 
and n max « 2 x 10 11 cm~ 2 . The second criterion gives 
the dashed curve in Fig. [T3l which is, on average, almost 
one order of magnitude larger than without a magnetic 
field. Note that the effect of the magnetic field is much 
stronger for large distances d. 

B. Qubit interaction at the interface: exchange 

One of the problems for quantum computation in 
doped Si arises from the lack of control of the exact po- 
sition of the donors. The main consequence of this is the 
indetermination of the value of the exchange between two 
neighboring donor electrons due to the theoretically pre- 
dicted oscillations of exchange with R, caused by valley 
interference effectsi&i^ One straightforward way to alle- 
viate this problem is to perform these operations at the 
interface^ where, as discussed in Sec. Ill B[ this degener- 
acy is lowered. Additionally, it would be much easier to 
control the qubit operations when the electrons are at the 
interface, similar to the successful experiments on dou- 
ble quantum dots in GaAs 2 ^ 2 ^^ and Si^ Note that 
the potential created by donor pairs (inset in Fig. [TS)) re- 
sembles very much a double quantum dot with the clear 
advantage that, in this case, the potential is produced 
exclusively by the Coulomb attraction of the donors and 
its exact form is known: V = 2{Q — l)/\/p 2 + d 2 . 

It is therefore of interest to determine the exchange 
coupling between donor electrons at the interface. As 




Inter-donor distance R(nm) 

FIG. 14: (Color online) (a) Overlap S between electron wave- 
functions at neighboring wells at the interface for three differ- 
ent values of d. (b) Exchange J calculated within the Heitler- 
London approximation. Values for S and J given here satisfy 



a first approach, we perform these calculations within 
the Heitler-London method. The validity and limita- 
tions of this method to calculate exchange in semicon- 
ductor nanostructures has been previously discussed by 
the authors^ The expression for the exchange within 
this approximation is 

25 2 e 2 
f - A 4 e S iri2 

- ^^(^(1)^.(2)12^(1) + -^H$£(2)*h(1)> 
1 S e S iri2 

(24) 

where $ Ljfl = g(x T Rf2,v), S = exp[-/3 2 (i?/2) 2 ] is the 
overlap, and e 2 /es\fi2 is the electron-electron interaction 
with r\2 the distance between electron (1) and electron 
(2). The first term in Eq. |24|) is the direct term and the 
second is the exchange term. 

In Fig. [TJ] we show the exchange J and the overlap 
S versus the inter-donor distance R for three different 
values of d. Note that the overall dependence of these 
two quantities with R is very similar, indicating that the 
behavior of J is closely related to the overlap^ 2 - J and 
S values are shown only for distances R > Rmm, defined 
in Sec. MI Al For a wide range of i?'s, J is of the same 
order as in GaAs double quantum dots where J's as low 
as ~ 10 neV have been measured^ 

Exchange control can be performed by applying a mag- 
netic field perpendicular to the interface, which reduces 
the wave-function radius, therefore decreasing the over- 
lap. The effect of the magnetic field on the exchange 
is well known)^ J decreases and eventually changes 
sign when the triplet is favored becoming the ground 
state. This is illustrated in Fig. fJ5] for d = 6.3 nm and 
R = 28 nm. At very large fields the singlet and triplet 
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Magnetic field (T) 

FIG. 15: (Color online) Modulation of overlap and exchange 
by a magnetic field perpendicular to the interface in the par- 
ticular case d = 6.3 nm and R — 28 nm. The double well 
potential for these particular values of d and R is depicted in 
the inset of Fig. [13] 



states become degenerate, as expected. Note that, in this 
case, the qualitative behavior of S and J with magnetic 
field is very different. 

IV. SUMMARY AND CONCLUSIONS 

Quantum computer architectures based on semicon- 
ductor nanostructure qubits have the key potential ad- 
vantage of scalability, which has led to the great deal 
of current interest in Si- and GaAs-based quantum com- 
putation. Silicon based spin qubits have the important 
additional advantage of extremely long spin coherence 
times (T 2 Sl ~ miliseconds or more) since isotopic purifi- 
cation (eliminating 29 Si nuclei) could considerable sup- 
press spectral diffusion induced electron spin decoher- 
encc" M£&L leading to T 2 Si - 100 ms, whereas electron 
spin coherence time is constrained to be rather short 
in GaAs quantum dot structures, XIp 3 ^ ~ 1 — 10 fis, 
since neither Ga nor As nuclei have zero nuclear an- 
gular momentum isotopes. However, the experimental 
progress in Si spin qubits has been very slow whereas 
there has been impressive recent experimental progress in 
the GaAs quantum dot spin qubitsi22i22i22i22i2I The main 
experimental advantage of GaAs quantum dot system has 
been the ease in the 1-qubit and 2-qubit manipulation be- 
cause the electrons near the surface can be effectively con- 
trolled by surface gates. By contrast, Si:P qubits are in 
the bulk, severely hindering experimental progress since 
electron manipulation in the bulk has turned out to be a 
difficult task. 



In this paper we show through detailed quantitative 
theoretical work how to control and manipulate qubits 
(i.e. both single electrons and two-electron exchange cou- 
pling) in a doped-Si quantum computing architecture^ by 
applying external electric and magnetic fields. In partic- 
ular, we have analyzed three main issues: (i) the times in- 
volved in the donor electron 'shuttling' between the donor 
and the interface of Si with (typically) SiC"2 have been 
found to be a few orders of magnitude shorter than the 
spin coherence times in Si, as required to allow for the 
necessary 'logic operations' and 'error correction' to take 
place; (ii) the existence of a well defined interface state 
where the electron remains bound and localized, so that it 
does not spread and form a 2DEG. This condition, which 
guarantees that electrons actually involved in a particular 
operation be taken back from the interface to donor sites, 
leads to a lower bound for the interdonor spacing, and 
consequently a maximum donor planar density; and (iii) 
the possibility of performing the two-qubit exchange gate 
operations at the interface, instead of around the donor 
sites as originally proposed^ Our results show that suffi- 
ciently large values of exchange coupling (~ 10 meV) 
can be achieved. 

Interface operations have several potential advantages 
over bulk operations, the most obvious one being that the 
read-out procedure would be simplified. A well known 
limitation of exchange gates for donor electrons in the 
bulk is the oscillatory behavior of the exchange coupling, 
which is due to the strong pinning of the six conduction 
band Bloch-function phases at each donor site, where the 
Coulomb potential is infinitely attractive^ This condi- 
tion is alleviated at the interface in two ways: first, the 
six valley degeneracy is partially lifted and, second, al- 
though the electrons at the interface remain bound to the 
donors, the binding potential is not singular; it is actu- 
ally equivalent to a quantum dot potential. Experiments 
on charge-qubit control in a double quantum dot at Si 
surface 20 indicate that the exchange oscillatory behavior 
may not be a severe problem for donor-bound electrons 
manipulated at the Si/Si02 interface. 

Our proposal combines the advantages of Si spin qubits 
(i.e. long T 2 time) with the structural advantages of 
GaAs qubit control and manipulation. We believe that 
the specific experiments we propose (and analyze in 
quantitative details) in this paper, if carried out, will go 
a long way in establishing the feasibility of a Si quantum 
computer. 
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